Can we measure the neutrino mass hierarchy 
in the sky? 



Raul Jimenez 1,2 , Thomas Kitching 3 , Carlos Pena-Garay 4 
and Licia Verde 1,2 

1 ICREA & ICCUB-IEEC, University of Barcelona, Barcelona 
08028, Spam. 

2 Institute for the Physics and Mathematics of the Universe 
(IPMU), the University of Tokyo, Kashiwa, Chiba, 211-8568, Japan. 

3 Institute for Astronomy, University of Edinburgh, Blackford Hill, 
Edinburgh EH9-3HJ, UK. 

4 Instituto de Fisica Corpuscular (CSIC-UVEG), Valencia, Spain. 

Abstract. Cosmological probes are steadily reducing the total neutrino mass 
window, resulting in constraints on the neutrino-mass degeneracy as the most 
significant outcome. In this work we explore the discovery potential of 
cosmological probes to constrain the neutrino hierarchy, and point out some 
subtleties that could yield spurious claims of detection. This has an important 
implication for next generation of double beta decay experiments, that will be 
able to achieve a positive signal in the case of degenerate or inverted hierarchy 
of Majorana neutrinos. We find that cosmological experiments that nearly cover 
the whole sky could in principle distinguish the neutrino hierarchy by yielding 
'substantial' evidence for one scenario over the another, via precise measurements 
of the shape of the matter power spectrum from large scale structure and weak 
gravitational lensing. 



1. Introduction 

In the past decade, there has been great progress in neutrino physics. It has 
been shown that atmospheric neutrinos exhibit a large up-down asymmetry in the 
SuperKamiokande (SK) experiment. This was the first significant evidence for a finite 
neutrino mass [T] and hence the incompleteness of the Standard Model of particle 
physics. Accelerator experiments [2j [3] have confirmed this evidence and improved the 
determination of the neutrino mass splitting required to explain the observations. The 
Sudbury Neutrino Observatory (SNO) experiment has shown that the solar neutrinos 
change their flavors from the electron type to other active types (muon and tau 
neutrinos) |1] . Finally, the KamLAND reactor anti-neutrino oscillation experiments 
reported a significant deficit in reactor anti-neutrino flux over approximately 180 km 
of propagation 5 . Combining results from the pioneering Homestake experiment [6] 
and Gallium-based experiments [7], the decades-long solar neutrino problem [8] has 
been solved by the electron neutrinos produced at Sun's core propagating adiabatically 
to a heavier mass eigenstate due to the matter effect [5]. 

As a summary, two hierarchical neutrino mass splittings and two large mixing 
angles have been measured, while only a bound on a third mixing angle has been 
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established. Furthermore the standard model has three neutrinos and the motivation 
for considering deviations from the standard model in the form of extra neutrino 
species has now disappeared [TTJ [T^] . 

New neutrino experiments aim to determine the remaining parameters of the 
neutrino mass matrix and the nature of the neutrino mass. Meanwhile, relic neutrinos 
produced in the early universe are hardly detectable by weak interactions but new 
cosmological probes offer the opportunity to detect relic neutrinos and determine 
neutrino mass parameters. 

It is very relevant that the maximal mixing of the solar mixing angle is excluded 
at a high significance. The exclusion of the maximal mixing by SNO [4] has an 
important impact on a deep question in neutrino physics: "are neutrinos their own 
anti-particle?". If the answer is yes, then neutrinos are Majorana fermions; if not, 
they are Dirac. If neutrinos and anti-neutrinos are identical, there could have been 
a process in the Early Universe that affected the balance between particles and anti- 
particles, leading to the matter anti- matter asymmetry we need to exist [13]. This 
question can, in principle, be resolved if neutrinoless double beta decay is observed. 
Because such a phenomenon will violate the lepton number by two units, it cannot 
be caused if the neutrino is different from the anti-neutrino (see |10) and references 
therein). Many experimental proposals exist that will increase the sensitivity to such 
a phenomenon dramatically over the next ten years (e.g., |14) and references therein). 

The crucial question we want to address is if a negative result from such 
experiments can lead to a definitive statement about the nature of neutrinos. Within 
three generations of neutrinos, and given all neutrino oscillation data, there are three 
possible mass spectra: a) degenerate, with mass splitting smaller than the neutrino 
masses, and two non-degenerate cases, b) normal hierarchy, with the larger mass 
splitting between the two more massive neutrinos and c) inverted hierarchy, with the 
smaller spitting between the two higher mass neutrinos. For the inverted hierarchy, 
a lower bound can be derived on the effective neutrino mass [10]. The bound for the 
degenerate spectrum is stronger than for inverted hierarchy. Unfortunately, for the 
normal hierarchy, one cannot obtain a similar rigorous lower limit. 

Neutrino oscillation data have measured the neutrino squared mass differences, 
which are hierarchical. Given the smallness of neutrino masses and the hierarchy in 
mass splittings, we can characterize the impact of neutrino masses on cosmological 
observables and in particular on the 

the matter power spectrum by two parameters: the total mass £ and the ratio of the 
largest mass splitting to the total mass, A. As we will show, one can safely neglect 
the impact of the solar mass splitting in cosmology. In this approach, two masses 
characterize the neutrino mass spectrum, the lightest one, to, and the heaviest one, 
M. 

Neutrino oscillation data are unable to resolve whether the mass spectrum consists 
in two light states with mass to and a heavy one with mass M, named normal hierarchy 
(NH) or two heavy states with mass M and a light one with mass to, named inverted 
hierarchy (IH). Near future neutrino oscillation data may resolve the neutrino mass 
hierarchy if one of the still unknown parameters that relates flavor with mass states 
is not too small. On the contrary, if that mixing angle is too small, oscillation data 
may be unable to solve this issue. Analogously, a total neutrino mass determination 
from cosmology will be able to determine the hierarchy only if the underlying model 
is normal hierarchy and E < 0.1 eV (see e.g., Fig[lJ. If neutrinos exist in either an 
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Figure 1. Left: constraints from neutrino oscillations and from cosmology in the 
m-T, plane. Right: constraints from neutrino oscillations (shaded regions) and 
from cosmology in the E-A plane. In this parameterization the sign of A specifies 
the hierarchy. 



inverted hierarchy or are denegerate, (and if the neutrinoless double beta decay signal 
is not seen within the bounds determined by neutrino oscillation data) , then the three 
light neutrino mass eigenstates (only) will be found to be Dirac particles. 

In this paper, we investigate whether cosmological data may positively establish 
the degenerate spectrum from the inverted hierarchy (or vice versa) . Our approach is 
to take cosmic variance limited surveys, rather than specifically planned experiments, 
so that we can determine if (even in the ideal case) cosmology can make any impact 
on this question. 

2. Massive Neutrinos and the Power Spectrum 

Massive neutrinos affect cosmological observations in a variety of different ways. For 
example, cosmic microwave background (CMB) data alone constrain the total neutrino 
mass E < 1.3 eV at the 95% confidence level [15]. Neutrinos with mass < leV become 
non-relativistic after the epoch of recombination probed by the CMB, thus massive 
neutrinos alter matter-radiation equality for a fixed Q m h 2 . After neutrinos become 
non-relativistic, their free streaming damps the small-scale power and modifies the 
shape of the matter power spectrum below the free-streaming length. Combining 
large-scale structure and CMB data, at present the sum of masses is constrained to 
be E < 0.3 eV [TB]- Forthcoming large-scale structure data promise to determine 
the small-scale (0.1 < k < 1 h/Mpc) matter power spectrum exquisitely well and to 
yield errors on E well below 0.1 eV (e.g., jTTl HHJ ) • Here, we assume the standard 
ACDM model and explore the changes in the matter power spectra due to the neutrino 
properties (mass and hierarchy). 

The effect of neutrino mass on the CMB is related to the physical density of 
neutrinos, and therefore the mass difference between eigenstates can be neglected. 
However individual neutrino masses can have an effect on the large-scale shape of the 
matter power spectrum. In fact, neutrinos of different masses have different transition 
redshifts from relativistic to non-relativistic behavior, and their individual masses 
and their mass splitting change the details of the radiation-domination to matter- 
domination regime. As a result the detailed shape of the matter power spectrum on 
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scales k ~ 0.01 h/Mpc is affected. In principle therefore a precise measurement of the 
matter power spectrum shape can give information on both the sum of the masses and 
individual masses (and thus the hierarchy), even if the second effect is much smaller 
than the first. 

We dehne the relation between the neutrino masses m and M and the parameters 
£ and A as 

NH : E = 2m + M A = (M-m)/£ (1) 
IH : £ = m + 2M A = (m - M)/£ (2) 

(recall that m denotes the lightest neutrino mass and M the heaviest). 

In Fig [1] we show the current constraints on neutrino mass properties in the m-A 
and £-A planes. While many different parameterizations have been proposed in the 
literature to account for neutrino mass splitting in a cosmological context P2B [5D1 [3T] 
here we advocate using the A parameterization for the following reasons. A changes 
continuously through normal, degenerate and inverted hierarchies; A is positive for 
NH and negative for IH. Finally, as we will show, cosmological data are sensitive to A 
in an easily understood way through the largest mass splitting (i.e., the absolute value 
of A), while the direction of the splitting (the sign of A) introduces a sub-dominant 
correction to the main effect. This parameterisation is strictly only applicable for 
£ > 0, but oscillations experiments already set £ > M > 0.05eV. 

It is important to note that not the entire parameter space in the £-A plane 
(or of any other parameterization of the hierarchy used in the literature) is allowed 
by particle physics constraints and should be explored: only the regions around the 
normal and inverted hierarchies allowed by neutrino oscillation experiments should be 
considered (see Fig[T]). 

To gain a physical intuition on the effect of neutrino properties on cosmological 
observables, such as the shape of the matter power spectrum, it is useful to adopt 
the following analytical approximation, as described in Ref. |20j . The matter power 
spectrum can be written as: 

where A^ is the primordial amplitude of the fluctuations, n s is the primordial power 
spectrum spectral slope, T(k) denotes the matter transfer function and D u (k,z) is 
the scale-dependent linear growth function, which encloses the dependence of P(k) on 
non-relativistic neutrino species. 

Each of the three neutrinos contributes to the neutrino mass fraction f v ^ where 
i runs from 1 to 3, 

'*-fe-™(a&v)Gs&) < 4 » 

and has a free-streaming scale kf s ^, 



Analogously, one can define the corresponding quantities for the combined effect of all 
species, by using £ instead of m Vi . Since we will only distinguish between a light and 
a heavy eigenstate we will have e.g., f u ,nn fv,Y,, fcfs,m, &fs,£ etc., where in the expression 
for fu,m one should use the mass of the eigenstate (which is the mass of the individual 
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neutrino, or twice as much depending on the hierarchy) while in k[ s m one should use 
the mass of the individual neutrino. 

The dependence of P(k) on non-relativistic neutrino species is in D„(h, z), given 

by 

D Vi (k,z)<x(l-f Vi )D(z) 1 -» (6) 

where k 3> kf B ,i {z) and pi — (5 — y25 — 24/^) /4. The standard linear growth function 
-D(z) fitting formula is taken from |22) . 

In summary there arc three qualitatively different regimes in fc-space that are 
introduced by the neutrino mass splitting 

D v (k,z) = D(k,z) fc<fc fs , m (7) 

D u (k,z) = (l-/v,m)I>(2) (1 ~ Pm) fcfs,m<fc< fc&.E (8) 

z) = (1 - /„, s )I^) (1 ~ Ps) fc > fc fs , s , (9) 

where the subscript m refers to the light neutrino eigenstate and S to the sum of all 
masses. 

This description is, however, incomplete: the transitions between the three 
regimes is done sharply in k while in reality the change is very smooth. In addition, the 
individual masses change the details of the matter-radiation transition which (keeping 
all other parameters fixed) adds an additional effect at scales k > fcf s ,£. 

In order to explore what constraints can be placed on A and £ for a given survey 
set-up we can use a Fisher matrix approach. The elements of F, the Fisher information 
matrix |23j . are given by 

where 9 and 7 denote cosmological parameters (and the Fisher matrix element's 
indices) and L denotes the likelihood of the data given the model. Marginalised errors 
on a parameter are computed as cr 2 (9) = (F~ 1 )gg . We can also calculate expected 
Bayesian evidence for cosmological parameters using the approach of Ref. |33] • In 
the case that we are considering we use the formula from [33] for the expectation value 
of the evidence, in this case the expected Bayes factor is simply the log of ratio of the 
Fisher determinants. 

Following Ref. [25] the Fisher matrix for the galaxy power spectrum is 

y ri /-fcmax dbxP(k,u)d]nP(k,u) 
F 6l = - 2 J ^ d, edUN ^ £ (11) 

with N = [nP(k,fi)/(nP(k,fi) + l)] 2 and V s is the volume of the survey. The 
integration over the projected angle along the light of sight Q /i is analytical and the 
maximum and minimum wavenumbers allowed depend on the survey characteristics 
with the constraint that A: max must be in the linear regime. The derivatives are 
computed at the fiducial model chosen. Throughout this paper we assume a fiducial 
model given by basic parameters of the standard LCDM cosmology [15] and the fiducial 
values for £ and A are then further specified in each case. 

Despite its limitations, the analytic description of the neutrino effect described 
above is extremely useful when performing an order-of magnitude calculation of an 
effect. Its corresponding Fisher matrix-approach forecasted errors indicate that while 

X As it is customary, fi denotes the cosine of the angle with respect to the line of sight. 
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Figure 2. Dependence of P(k) on the parameter A at z = 0, for fixed £ and 
several values of A. The dependence is expressed as fractional variation in P(k) 
for a unit variation in A. For this value of the total mass S, normal (inverted) 
hierarchy correspond to A ~ 0.05 (A = —0.05). 

E can be constrained tightly, nearly-ideal, full-sky, cosmic- variance dominated surveys 
will be needed to obtain promising errors on A. However, we find that the analytical 
approximation above overestimates the neutrino effects on the P(k) and therefore 
under-estimates the forecasted errors, by factors of ~ few in the regime of interest 
(E < 0.3 eV, A along the NI and IH lines). In what follows we therefore use the 
publicly available CAMB code 26 to compute the matter power spectrum. 

In Fig. [2] we show the dependence of P(k) on the parameter A at z = for fixed E 
and hxed cosmological parameters. The dependence is shown as the fractional change 
of the matter power spectrum for a unit change of the parameter A. This quantity is 
then fed directly to the Fisher matrix (see Eq. Hip . In order to compute reliably the 
above derivatives, CAMB needs to be run at the highest precision settings, with fine k 
sampling and taking care that interpolations procedures in-built in the code do not 
introduce a spurious signal. 

For the value of the total mass E = 0.25 eV, adopted in Fig. [2 normal (inverted) 
hierarchy correspond to A ~ 0.05 (A ~ —0.05), indicating that the effect of the 
neutrino mass splitting on the P(k) is at the ~ 0.2% level. The dependence of P(k) 
on A at k > 0.1 h/Mpc arises because even for a fixed E the individual masses affect 
the tail of the energy distribution of the relativistic species and thus matter-radiation 
equality. Note that <91nP/<9A changes sign with A and there is a location, A = 0, the 
degenerate case, where P(k) shows no dependence on A. 

To understand the meaning and implications of this let us consider that the error 
on A is directly proportional to A^ 2 = — 2(lnL — lnigduciai) where ^fiducial denotes 
the Likelihood for the fiducial model. In addition we can write 

A x 2 cx / k 2 [P(k,A)-P(k,A Muclal )] 2 dk. (12) 

•/fcmto 
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Figure 3. A% 2 as a function of the degeneracy parameter A for a fixed total 
neutrino mass S (and fixed cosmology). This is a section along a £=constant line 
of Fig. [T]of the quantity — 21nL as it would be seen by a Fisher matrix approach 
for a IH fiducial model. The vertical lines show the location of the normal and 
inverted hierarchy. Note the bimodal distribution of the In L surface, which makes 
the determination of the hierarchy from measurements of the shape of the power 
spectrum extremely challenging. The Ax 2 normalization matches that achievable 
from an ideal weak lensing survey as described in the text. 

This quantity is shown in Fig. [3l where the normalization has been chosen to 
match the constraints achievable from an ideal full sky weak lensing survey as the one 
considered in Sj3] 

Fig. |3] shows that the likelihood surface is bimodal: for example, for a fixed 
cosmology and fixed E, if the fiducial model is the inverted hierarchy A < 0, there 
is a corresponding value of A > (normal hierarchy), consistent with the neutrino 
oscillations in the allowed region, with a P(k) virtually indistinguishable from the 
fiducial model. 

The bimodality of the likelihood surface also implies that the Fisher matrix 
approach to forecasting errors need to be applied with care before it can interpreted 
in terms of distinguishing the hierarchy. The curvature of the likelihood around the 
fiducial model gives the formal error on A, and this error may be much smaller 
than the distance between IH and NH A values. But this could be interpreted as 
a determination of the hierarchy if and only if the likelihood had a unique maximum, 
which is not the case here. This subtlety has not be noticed in the literature 
before where Fisher error-estimates for parameters describing neutrino hierarchy were 
presented. In general, errors have been computed around one or more fiducial 
models and were sometimes found to be smaller than the distance between normal 
and inverted. We point out here that this cannot be directly interpreted as being 
sufficient to distinguish the hierarchy (it is a necessary but not sufficient condition if 
the likelihood is multi-peaked). 

A more detailed inspection of Fig. [3] also indicates that the A^ 2 between the two 
minima (maxima of InL) is not exactly zero, but it is very small, and that the location 
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Figure 4. LSS (left) and Weak Lensing (right) forecasts for neutrino mass 
parameters E and A. We assume the LSS survey consists of a comoving volume 
of 600 Gpc 3 at z = 2 and 2000 Gpc 3 at z = 5. The Weak Lensing survey covers 
40,000 sq. deg. with a median redshift of 3.0 and a number density of 150 galaxies 
per square arcminute. Several fiducial models (E,A) were used to derive by Fisher 
matrix approach the expected errors. The upper panel show the l-a errors on 
A and E, the errors in E are so small that are barely visible. The lower panel 
shows the expected evidence ratio between the normal and inverted constraints 
as a function of neutrino mass. The dashed line shows the InB = 1 level: in 
Jeffrey's scale lnS < 1 is 'inconclusive' evidence, and 1 < InB < 2.5 corresponds 
to 'substantial' evidence. 

of the second minimum (assuming a fiducial IH) does not coincide with the central 
value of the oscillations-experiments regions. The evidence ratio can the be used to 
quantify wether a given survey set up could distinguish the two cases. 

The philosophy of the rest of the paper is therefore: "can cosmology in the cosmic- 
variance limit, and for an ideal experiment, distinguish the neutrino heirarchy?" or 
in other words, "is there enough information in the sky to measure the neutrino 
hierarchy? " 

3. Forecasted Constraints from Large Scale Structure 

Here we explore what constraints can be placed on A and £ from ideal, cosmic 
variance-dominated future surveys probing the shape of the matter power spectrum. 
The two probes of large-scale structure (LSS) we consider are the matter power 
spectrum itself and weak lensing. 

We also compute the Fisher matrix of a CMB experiment like Planck [j| in 
order to help break degeneracies in the cosmological parameters when determined 
only by the power spectrum, or weak lensing. Therefore our final Fisher matrix is 
F = Fp(k),WL + Fcmb- We compute the combined Fisher matrix for variations in 
the following cosmological parameters: n s , a s , ft v h 2 , A, Z, ftbh 2 , fi c /i 2 , h, A s , where a s 
denotes the running of the power spectrum spectral slope and Z is related to the 
optical depth to the last scattering surface via Z = cxp(— 2r); Sl„ is related to £ via 
£ = 9451^, (eV). The reported errors on £ and A are marginalized over the other 
cosmological parameters. The marginal errors for A and £ are shown in Fig. H] the 

§ www.sciops.csa.int/PLANCK/ 
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left panel is for a direct P(k) measurement approach and the right panel is the weak 
lensing approach. 

Because we are interested in answering the question: "is there enough information 
in the sky to measure the neutrino hierarchy?" we have chosen survey parameters that 
are ambitous cosmic variance-limited surveys. For the parameter points shown in the 
left panel of Fig. [4] we have assumed a survey that covers the full sky 40,000 square 
degrees and maps the positions of galaxies up to z — 2 corresponding to about 600Gpc 3 
comoving volume and maps the 21cm-HI up to z = 5, corresponding to about 2000 
Gpc 3 comoving volume. We also assumed a high number density of galaxies so that 
we work in the cosmic variance-dominated regime (nP ^> 1). Galaxies are expected 
to be a biased tracer of the dark matter distribution, here we assume the bias to be 
scale and redshift-independent and thus not to affect the recovery of the shape of the 
matter power spectrum. 

HI surveys [27] target the hyperfine transition in the hydrogen atom, which in 
the rest-frame emits a photon in the radio wavelenghts (21 cm). Therefore they 
survey the amount of neutral hydrogen in the universe. Because most galaxies and 
dark matter overdensities contain neutral hydrogen, such surveys provide the most 
un-biased indirect tracer of the dark matter distribution in the Universe. Further, in 
this frequency band, the radio spectrum is featureless with the only line being the 21 
cm one, its observed frequency yielding a redshift and thus the radial distance of the 
emitter. Thus, an imaging survey automatically gives a three dimensional map of the 
HI distribution. The main challenge facing the HI surveys is the contamination by 
foregrounds [27]. For the characteristics of the survey we have followed the numbers 
given in [28] which yield to bias of 1 and negligible shot noise. 

The survey considered is certainly a challenging one, but our calculations indicate 
that a cosmic variance-limited galaxy and HI survey can provide enough information 
to determine the neutrino hierarchy. We find that such a survey could constrain the 
total sum of neutrino mass with extreme accuracy O(10~ 5 ). We also find that if the 
total neutrino mass is smaller than 0.15 eV, then the IH could be distinguished from 
the NH through an evidence criteria centered on each peak in A. 

Weak lensing is the effect where the path of photons propogating from a galaxy 
are distorted by intervening mass concentrations. The amount of distortion depends 
on the density and distribution of the mass. For an individual galaxy image the weak 
lensing effect is to induce a change in ellipticity or 'shear'. By using redshift and 
shear measurements from every galaxy, information on the growth of structure and 
the geometry of the Universe can be extracted from 3D cosmic shear observations. 
Here we will use the 3D cosmic shear approach [29] [30] [31] where the full 3D shear 
field is characterised using 3D spherical harmonics and the Fisher matrix methodology 
of [31j . In line with the cosmic- variance limited approach of this article, we assume 
a large, cosmic-variance limited weak lensing survey covering 40,000 square degrees, 
to a median redshift of 3 with 50 galaxies per square arcminute. On the right panel 
of Fig. [H we show the marginalised constraints A and E, for this cosmic- variance 
limited survey. We find that the sum of neutrino mass is constrained to extreme 
accuracy 0(1O -6 ). As the neutrino mass decreases the constraints on the IH and NH 
become smaller and for massess below ~ 0.15 eV the evidence ratio for the IH and NH 
constraints (lower panels of Fig. HJ would become substantial (in a Jeffrey's scale), 
allowing for the neutrino heirachy to be distinguished. This is again a very challenging 
survey, and acts to show highlight how demanding the measurement of neutrino mass- 
splitting can be; however by using shear measurements from Euclid [32] or LSST [17] 
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we may hope to approach this regime. 

The degeneracies between E and A are small, and the very small constraint on 
E results in the constraints being effectively un-correlated in the E-A plane. We note 
that the constraints on A around the IH and NH peaks are tighter for weak lensing 
than LSS, this is due to lensing providing constraints on both the geometry and the 
growth of structure, which provides a smaller raw constraint and a more orthogonal 
constraint to the CMB resulting in smaller errors. Interestingly, even though the 
weak lensing constraints on A are smaller than for the power spectrum, the evidence 
ratio is comparable, because, due to the multi-dimensional degeneracy directions, a 
naive correspondence between error-bars and evidence is not applicable (it is to a first 
approximation the difference between the two error bars that is important). 

Note that the evidence calculation explicitly assumes two isolated peaks, and so is 
only applicable when the fiducial points are seperated by multiple-sigma. As a result 
of this, the evidence calculations may be slightly optimisic for large masses. However, 
for S < 0.2 eV, the x 2 difference between the two minima becomes noticeable as well 
as the shift between the location of one of the two minima and the central A value 
for the oscillations experiments (which induces an additional \ 2 difference). While 
this information is not fully accounted for in a Bayesian approach to forecasting the 
evidence, it may be included at the moment of analyzing the data, using different 
approaches such as the likelihood ratio, and may slightly improve the significance for 
the hierarchy determination. 

While we have used the oscillation results to center the Fisher and evidence 
calculations on the NH and IH, combining the oscillation experiments constraints will 
not improve the evidence; in fact, oscillation experiments give symmetric errors on A 
(i.e. they do not depend on the sign of A). 

4. Conclusions 

The shape of the matter power spectrum contains information, in order of decreasing 
sensitivity, about the sum of neutrino masses, the amplitude of the mass splitting 
and the hierarchy (i.e., the mass splitting order). We have introduced a novel 
parameterization of the neutrino mass hierarchy, A, that has the advantage of changing 
continuously between normal, degenerate and inverted hierarchies and whose sign 
changes between normal and inverted. The absolute value of A describes the maximum 
mass difference between the eigenstates. We stress that, current constraints from 
neutrino oscillations have ruled out large part of the parameter space given by the 
sum of the masses and the A parameter, leaving two narrow regions: for a fixed value 
of the total mass, the value of A for the normal hierarchy is related to that of the 
inverted one and Anh — — Ajh (but, in detail, Anh ^ |A/#|). It is the allowed 
region that cosmology should explore. 

We found that the information about A accessible from the power spectrum 
shape yields a degeneracy: parameters values A and —A yield nearly identical power 
spectra and therefore that the likelihood surface in A is bimodal. This was not noted 
in the literature before and not taking this into account when using the Fisher matrix- 
approach to forecast future surveys performance may lead to spurious indications of 
a surveys ability to determine the hierarchy. 

Detecting the signature of the hierarchy in the sky is therefore extremely 
challenging, and therefore we asked: "can cosmology in the cosmic-variance limit, 
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Figure 5. Role of cosmology in determining the nature of neutrino mass. Future 
neutrinoless double beta decay (Oufibeta) experiments and future cosmological 
surveys will be highly complementary in addressing the question of whether 
neutrinos are Dirac or Majorana particles. Next generation means near future 
experiments whose goal is to reach a sensitivity to the neutrinoless double beta 
decay effective mass of 0.01 eV. We can still find two small windows where this 
combination of experiments will not be able to give a definite answer, but this 
region is much reduced by combining Qufibeta and cosmological observations. 



and for an ideal experiment, distinguish the neutrino heirarchy?" or in other words, 
"is there enough information in the sky to measure the neutrino hierarchy?" To address 
these questions we have considered ideal, full-sky, cosmic variance-limited surveys and 
found that substantial Bayesian evidence (lnB > 1) can be achieved. Are such a 
surveys feasible in the next 5-10 years? There are two candidates for such surveys : a 
full extragalactic survey in the optical/infrared like Euclicfj] [32 and a full 21cm survey 
by the SKAIjI. Each of these surveys is scheduled to start operations by 2018. Euclid 
will make an all sky Hubble-quality map for weak lensing and will directly trace the 
dark matter using this technique; whilst the cosmic variance limited survey we consider 
here is ambitous with respect to this survey these result serve as a qualitative measure 
of this surveys expected performance (costraints should be only a factor < 1.5 larger 
at worst). Euclid will also target emission line galaxies up to z ~ 3 (therefore these 
galaxies will have bias parameter close to 1) however nP, quantifying the the ratio 
of the signal to shot noise, will be only slightly above 1. The 21cm surveys provide 

|| http://sci.esa.int/euclid 

If http://www.skatelescope.org 
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the most un-biased indirect tracer of the dark matter distribution in the Universe and 
have negligible shot noise. 

For the degenerate and inverted mass spectra, the next generation neutrinoless 
double beta decay experiments can determine if neutrinos are their own anti- 
particle. For the normal hierarchy, the effective electron-neutrino mass may even 
vanish. However, if the large-scale structure cosmological data, improved data on 
the tritium beta decay, or the long-baseline neutrino oscillation experiments establish 
the degenerate or inverted mass spectrum, the null result from such double-beta decay 
experiments will lead to a definitive result pointing to the Dirac nature of the neutrino 
mass. This is summarized in figure [5j 

If the small mixing in the neutrino mixing matrix is negligible, cosmology might 
be the most promising arena to help in this puzzle. Our work shows that depending 
on the total neutrino mass, there might be substantial evidence by cosmological data 
to infer the neutrino hierarchy. 
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